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Continuous-Time Quantum Monte Carlo Study 
of Local Non- Fermi Liquid State in the Multichannel Anderson Model 
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The impurity Green's function Gj in the local non-Fermi liquid state is evaluated by means of 
the continuous-time quantum Monte Carlo method extended to the multichannel Anderson model. For 
N — M (where N and M are numbers of spin components and channels, respectively), Gf is expressed 
as — ImG/(aj + iO) = c — | x ' 2 , and the zero-frequency value c depends only on N (= M). A corre- 
sponding impurity self-energy at low frequencies is composed of two parts: a resonance term related to c, 
and a non-Fermi liquid term proportional to l^] 1 ^ 2 . The characteristic energy scale is discussed in terms 
of the non-Fermi liquid term in the self-energy. 
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1. Introduction 

The multichannel Kondo effect is a typical example that 
leads to a local non-Fermi liquid ground state. 1 ' It has been 
recognized that the peculiar low-temperature behaviors ob- 
served in uranium compounds and metals with uranium im- 
purities are due to the two-channel Kondo effect. 2 ' This 
kind of non-Fermi liquid state has been investigated from a 
more general point of view based on models generalized to 
SU(iV) ® SU(M) symmetry. 3 ' Then, their critical nature has 
been discussed extensively. 4 ' 5 ' 

Regarding the (single-channel) Kondo problem, the Ander- 
son Hamiltonian gives clear insight: 6,7 ' the ground state is 
connected to that in the non-interacting limit. In this anal- 
ogy, the multichannel Kondo effect can be addressed based 
on an Anderson Hamiltonian. 8 ' The inclusion of the impurity 
charge degree of freedom enables us to describe the local dy- 
namics via the impurity Green's function. We thus consider 
the SU(7Y) <g> SU(M) multichannel Anderson model given 
by 2 ' 

ka.fi a 

+F^(X /1 ,_ Q c V +ri.c.). (1) 

CtfX 

The (pseudo-)spin index // and channel index a run over N 
and M components, respectively. The f 2 state forms a 
channel singlet (—a denotes the counterpart of a), and the f 1 
state | a) has the energy E cyi relative to The Hilbert space 
of / states is restricted to |a) and by using the X-operators 
X ln i — |7)(7'| with 7 = a, fi, on which ^ 7 X ln = 1 

is imposed. c QM = iV ' ^2 k Cfc QM with Nq being number 
of sites. The M -channel Coqblin-Schrieffer model is derived 
from the Hamiltonian (1) as a localized limit V 2 , E CK —> oo 
with V 2 /E cx fixed. Exact thermodynamics of the model (l) 9 ' 
as well as the localized limit 5 ' has been derived. 

Concerning the dynamical properties, a two-channel case, 
N = M = 2, has been clarified by the numerical 
renormalization group 10 ' and by an exact method. 11 ' Gen- 



eral cases have been investigated by perturbational treat- 
ments. 3 ' 12 ' 13 ' In this paper, we numerically investigate the dy- 
namical properties of the multichannel Anderson model. To 
this end, we develop an algorithm based on the recently de- 
veloped continuous-time quantum Monte Carlo (CT-QMC) 
method, 14 ~ 16 ' which is explained in the next section. We show 
numerical results for the impurity Green's function and self- 
energy in §3. 

2. CT-QMC for the multichannel Anderson model 

We study the model (1) by the CT-QMC, which evalu- 
ates a perturbation expansion stochastically. In the present 
case, we adopt the hybridization expansion. 15 ' Since the non- 
perturbative part is diagonal with respect to a and fi, the ef- 
ficient algorithm using a 'segment' picture is applicable by a 
slight modification. Figure 1 shows a diagram of a configu- 
ration of order V 6 . Spin states \i{ and channel states on ap- 
pear alternately, which are hereafter referred to as segment 
and anti-segment, respectively. In general, a configuration of 
order V 2k is represented by qk = {n, r/, a,, /Xj}. The trace 
over the local states is thus taken into account graphically. On 
the other hand, the trace over conduction electrons is evalu- 
ated based on Wick's theorem. A Monte Carlo sampling is 
performed in the configuration space composed of k and ■ 

We perform the following update processes: (i) addi- 
tion/removal of a segment or an anti-segment, and (ii) ex- 
change of spin or channel indices. Fig. 2(a) shows the ad- 
dition of a segment. The index /i of the segment is randomly 
chosen, and accordingly the update probability differs from 
that in ref. 15 by a factor of N. When either N or M is larger 
than 2, the ergodicity is not satisfied only by process (i). For 

I M-3 t I Mg t I Hi I 

P *3 T 3 X 2 X 2 T 'l X l 

Fig. 1. Diagrammatic representation of a configuration of order V 6 . The 
outgoing and incoming allows indicate creation and annihilation of con- 
duction electrons, respectively. 
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Fig. 2. Update processes: (a) addition of a segment, and (b) exchange of 
spin indices. 




Fig. 3. Examples of diagrams which cannot be reached without the 'ex- 
change' process. The spin and channel components are labeled as 1, 2, 
• • ■ , and a, b, ■ ■ ■ , respectively. 



example, configurations shown in Fig. 3 cannot be reached. 
This problem can be solved by introducing a process shown 
in Fig. 2(b), which exchanges the spin indices. We perform a 
similar update to exchange the channel indices as well. 

In the simulation, we observe negative weight configura- 
tions for N = M > 2. However, since their contribution is 
less than 10% in the parameter range shown in this paper, the 
sign problem has little effect on the simulation. 

3. Numerical Results 

In this paper, we restrict ourselves to N = M. We use a 
rectangular density of states p(e) = (1/2D)6(D — |e|) for 
conduction electrons with D = 1. We fix NV 2 = 0.12 so 
that the exponent of the Kondo temperature is the same for 
different N. The width of the localized state A = TrV 2 p(0) is 
A ~ 0.094 at most (for N = 2), and therefore the effect of 
finite band width may be neglected. 

3.1 Green's function 

We first show results for the single-particle Green's func- 
tion Gf, which is defined in the restricted Hilbert space by 



G f (ie n ) = - / dr(X^(T)X^ 
Jo 



(2) 



where e„ = (2n + 1)ttT is the fermionic Matsubara fre- 
quency. At high frequencies, Gf follows G/(ie„) ~ a/ie„ 
with a < 1, since the Hilbert space is restricted. The a varies 
between 1/N and 1/M depending on E cx , and in a special 
case of JV = M, a = l/N. 

In Fig. 4, — ImG/(ie„)A is plotted against y/e^ for E cx = 
0. For all N = M, Gf is expressed as — ImG/(ie„) = 
c — b^/e^ at low frequencies. Hence, G/(z) is non-analytic 
at z — > +i0, and the spectrum — ImG/(cj + iO) on real fre- 
quencies exhibits a cusp structure expressed by c — 6'lwl 1 / 2 , 
which has been reported for N = M — 2. 410,11) The value 
c at e„ — > +0 decreases with increasing N. From Fig. 4 and 
an analogy with the Friedel sum-rule in the Fermi liquid, we 
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Fig. 4. (Color online) The imaginary part of the Green's function Gf(ie n ) 
for N = M, NV 2 = 0.12, E cx = and T = 0.0001. The lines show 
sin 2 (7r/27V). 



conjecture the following relation: 

-ImG/(+i0) 



1 . 2 f n ^ 



. sin" I , 

A V27V/' 



(3) 



which is indicated in Fig. 4. Eq. (3) includes the result for 
N = M = 2, 4 ' 1(U1) 1/(2A), and reduces to the result in 
the non-crossing approximation, 3 ' 7r 2 /[(7V + M) 2 A], in the 
limit N = M 3> 2. For E cx = 0, the particle-hole symmetry 
leads to ReG/(ie„) = 0, meaning that the phase shift tf> of 
conduction electrons at cj = is fixed at = tt/2 irrespective 
of the value of N = M. Hence, the sine factor in eq. (3) is 
not connected with 4> but is due to the imaginary part of the 
self-energy. 

3.2 Self-Energy 

We discuss the self-energy in the restricted Hilbert space. 
The self-energy S/(ie n ) in the ordinary definition is given by 

Gf{ie n ) = -. — — — , (4) 

ie„ - ef - 1 (ie„J - hf{ie n ) 

where T(ie n ) = N^ 1 2~2k V 2 /(ie n — In the restricted 
Hilbert space, S/(ie„) diverges according to S/(ie n ) ~ 
ie„(l — 1/a) at e n — > oo, since G/(ie„) ~ a/ie n with a < 1. 
Although this divergence does not produce any problem with 
analysis of low-energy properties, it is not convenient in prac- 
tice. Thus, we define an alternative self-energy S f as follows: 

Gf(ie n ) = -. : ° (5) 

ie„ - ef - al (ie„) - L/(ie„J 

12 f is related to the ordinary self-energy E/ by S/ = aE / + 
ie„(l — a) — (if — aef), and converges in proportion to l/ie n 
at high frequencies. By using E f, for example, the renormal- 
ization factor z in the Fermi-liquid state is evaluated as 

Z=[l- 3ImE / (ie„)/ae„]7 n 1 ^ +0 
= a[l - 9ImE/(ie n )/9e n ]~ 1 _ > . +0 ee az. (6) 

z stands for a quasi-particle weight within the restricted 
Hilbert space, and accordingly E/ in eq. (5) may be a rea- 
sonable definition. 

Figure 5(a) shows — ImE/(ie ra ) divided by aA as a func- 
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Fig. 5. (Color online) (a) The imaginary part of the self-energy S^(ie„) 
defined by eq. (5). The lines show cot 2 (7r/2iV). (b) Sf(ie„) = 
T,f(ie„) — £ar(ie n ) on a log-log scale. The parameters are the same as 
in Fig. 4. 

tion of ^fe^. Similarly to G/(ie n ), E/(ie„) includes a term 
proportional to ^JT^ in the limit e„ — > +0, and converges to 
a finite value. To separate the zero-frequency value from 
we introduce a parameter £ as follows: 

^/( ie n) = £«r(ie„) + Sf(ie n ). (7) 

£ is determined so that ImS/(+iO) = 0. For N = M, noting 
that ReG/(ie„) = 0, we obtain from eq. (3) 

e=cot 2 (^). (8) 

In the case of N = M = 2, Sf(ie n ) eventually corresponds to 
the self-energy discussed in refs. 10 and 1 1 . Figure 5(b) shows 
Sf(ie n ) on a log-log scale. We can clearly see the power-law 
behavior — Im«S/(ie n ) oc lenl 1 ^ 2 , which means — lxaSf(uj + 
i0)cx|w| 1 / 2 . 

3.3 Effect of Level Splitting E cx : Energy Scale 

So far, we have examined E cx = 0. We now discuss the 
effect of £ cx . In refs. 10 and 1 1, it is reported for N = M = 2 
that ImG/(+i0) and ImS/(+i0) do not depend on E ex . We 
have confirmed for N — M > 2 that eqs. (3) and (8) hold up 
to E ex = 0.3 within numerical accuracy. The finite value of 
E ex causes an asymmetry of the cusp keeping the value at cj = 
0: c + x / 2 changes into c + [b + 9(uj) + 

As E cx increases, the energy scale becomes smaller. We 
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Fig. 6. (Color online) A characteristic energy scale Tq defined in eq. (9) 
as a function of E cx . To is evaluated from Sj(ieo) at T = 0.0005 
(open symbols) and T = 0.00025 (closed symbols). The lines show 
exp[-E SK /NV 2 p(0)]. 



define a characteristic energy scale To in terms of Sf by 

-lmS f (ie n )/aA - (e„/T ) 1/2 , (9) 

in the limit e n — > 0. Because To may be defined with an ar- 
bitrary factor, we shall discuss only its exponent. In Fig. 6, 
we show To as a function of E ex . Tq follows To oc Tk oc 
exp(-l/g) with g = NV 2 p{0)/E cx for E cx > 0.15, namely 
g < 0.4. We conclude that the exponent of the energy scale of 
the non-Fermi liquid self-energy agrees with the Kondo tem- 
perature Tk in the corresponding single-channel model. 

4. Summary 

We have presented the impurity Green's function G/(ie„) 
and the self-energy E/(ie„) in the non-Fermi liquid state 
using the CT-QMC extended to the multichannel Anderson 
model. For = M, G/ and £/ are non-analytic at oj = 
as |oj | 1 / 2 . The zero-frequency spectrum ImG/(+iO) does not 
depend on the excitation energy E cx , and correspondingly 
ImSj(+i0) has a finite value. These values depend only on 
N (= M), and seem to be expressed as eqs. (3), (7) and (8). 
An analysis of general N, M is left for future work. 
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